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A phenomenon analogous to the conical refraction well-known in the crystalooptics and crys- 
taloacoustics is considered for the magnetohydrodynamical waves in a coUisionless plasma with 
anisotropic thermal pressure. Imposing the most general condition for the existence of the phe- 
nomenon, angle of the conical refraction is calculated which appeared to be dependent on the ratio 
of the Alfven velocity and sound speed measured in the perpendicular direction in respect to the ex- 
ternal magnetic field. Feasible ways of experimental demonstration of the phenomenon are discussed 
and a novelty brought by the general consideration is outlined. 
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I. INTRODUCTION 



Recently, author of Ref.[l] showed, theoreticaUy, that in the magnetized plasma with anisotropic thermal pressure 
[viz., a medium where linear magnetohydrodynamic (MHD) waves are described by well-known double adiabatic theory 
^ of Chew, Goldberger and Low [2]] there is a possibihty of observation of the phenomenon of the conical refraction. 
T-H The similar kind of phenomenon was observed long time ago (see, for example hteraturc in Ref.[l]) for the hght and 
> acoustical waves in the biaxial crystals and therefore, it is very well established. For instance, in the biaxial crystals 
■ observation of the optical conical refraction is possible if the parallel, narrow beam of natural (non-polarized) light 
falls along the optical axis (binormal) of wave normals on a plate cut so that its two parallel faces are perpendicular to 
the binormal. While propagating through the crystal the beam will form an empty cone with permanently changing 
, linear polarization. Thus, producing a circle on the photographic plate placed after the crystal. 

In the Ref. [1] there were obtained the necessary conditions which ensure existence of the conical refraction in the 
medium under consideration. However, further investigation revealed that it is possible to impose much less restrictive 
conditions fulfillment of which ensures occurrence of the phenomenon. In the light of possible interest which plasma 
^ physics community may have to corroborate experimentally existence of the phenomenon, author find it worthwhile 
^ to outline here these novel conditions. 

In the section II we derive the basic equations and outline new, general conditions for the existence of the phe- 
O ^, nomenon. While in the section III we pay more emphasis to the physics of the effect; discuss possible ways of its 
experimental corroboration along with its feasible applications for plasma diagnostics; and finally we address a novelty 
brought by the general consideration employed in the present paper. 



, II. MAIN CONSIDERATION 

Let us consider the most general form of the dispersion relation for the MHD waves in magnetized plasma with 
anisotropic thermal pressure: 
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The latter equation can be derived from the dispersion relation obtained in Ref . [3] while considering its nonrelativistic 
limit (i.e. neglecting displacement current and ultrarelativistic temperature effects) , but retaining a freedom in choice 
of the polytropic indexes 6 and a. Note, that in the derivation of Eq.(l) there were used following equations of state: 



p± = const 



(2) 



P|| =const|,p(-| 



2a 



(3) 



We gather from the Eqs.(2) and (3) that when p|| = p± = p and a — 6 — the equations of state reduce to conventional 
one for isothermal, ideal gas P^p. For anisotropic pressure case, when py ^ p± and a = S = 1 they represent usual 
Chew, Goldberger and Low's equations of state [2]. Of course, taking into due account the two above limiting cases 
{a — S — Q and a — 6 — 1) in Eq.(l) results in the correct expressions for the dispersion relations for the linear 
MHD waves in isotropic and anisotropic pressure plasmas accordingly, which can be found in any textbook of plasma 
physics. 

For the sake of convenience we use further following notations [1]: 



Bo 



C_L = 4 / 




where a stands for the Alfven velocity, c± and C|| are the sound speeds measured in perpendicular and parallel 

directions in respect to the uniform, external magnetic field Bq. 
Now we can rewrite Eq.(l) in the following form: 



c^il + S sin^ 0) + 2acl cos^ 9 



±sj {a? + ci(l + ^ sin^ 6) - 2{a + l)c\ cos^ 6*)^ + Ac]_ sm^Ocos^ 6 



(4) 



where V = u/k denotes the phase velocity of the waves. 

In Ref.[l] after this step, in order get desired result (theoretical demonstration of the existence of the conical 
refraction) author required fulfillment of the following relation between C|| and c±: c± = \/2c|| . Furthermore, according 
to Ref.[l] final necessary condition for occurrence of the effect was obtained as: a = c_l = V^cy. However, now let us 
restrict ourselves with much more weak condition: 



+ ci - 2(a + 1)4 = 0, or cj 



2(a + l) 



(5) 



Substituting in Eq.(4) the latter expression for cf, and introducing x and if) as 



X = + (1 + S)cj_ and V = 



2 _ aa"^ - [d{a + 1) - a]c^ 



a + 1 



we can obtain: 



+ V'^ cos^ 9 ± + 2(2c| - x^) cos2 9 + {x'^ - 4ci) cos^ 9 



(6) 



Possessing the formula for the phase velocity of the wave it is easy to derive expression for the group velocity u via 
straightforward, but time consuming calculations: 

^ dw ->^^ av 

u = — = NV + k— = 
dk dk 



V 



[cos9/2V^] {^'^T{9) + (2ci - x'^) + (x^ - 4ci) cos^ 9) 
1 - C/V^ 



A+ 
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1 - C/V^ 



(7) 



where we introduced following notations: unit vectors A'' = k/\k\ and A = So/|-Bo|; T[6)=2V — (x + V' cos &) and 

C=(2a + l)[a2 + ci(l + (Jsin^ 0)] (^^^ cos^ ^ 



-(2a + 1) ( ° ^^^J cos^ e - cl sin^ ^ cos^ Q. 

Here we must note that in Eq.(7) we used biquadratic equation 

- + V'^ cos^ Q)V'^ + C = 0, 

solution of which dispersion relation (6) actually is. 

Let us describe a direction of the group velocity u by an angle between u and Bq (i.e., between u and ^4). For 
the tan we may write then 



tant^ ■ 



u X A 
u- A 



(2y2 _ ^2 ^,Qg2 _ (2c4 _ ;^4) ^Qg2 6» - (x^ - 44) cos" 6* 

(21/2 ^ ^2 sin2 ^g,^ _^ ^2c| - x'') sin^ + {x^ - 4ci) cos2 6» sin^ 



tan^. 



(8) 



Note, that T{0) when 6-^0. 

The latter equation determines the direction of the group velocity u whereas its modulus can be calculated by 
expression: 



V/ cos{ip - 0) 



(9) 



It is worth of remarking that the Eq.(9) is a consequence of the following expression u ■ N = V which can be readily 
shown by multiplying Eq.(7) scalarly by N. 

Now, it can be easily shown that the expressions for group velocity and taiup [Eqs. (7) and (8) correspondingly] 
exhibit uncertain behaviour as 0/0 when ^ Q. Removing of the uncertainty (for the sake of brevity we omit, here, 
several steps of this task) results in: 



tanwf+O) = limtan(p(0) = ± ( | ^ o 



= ± 



a + 1 \ p2 



2a + 1 



(10) 
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2a + 1\ J a'^+cl 
a + 1 



(11) 



In Eq.(lO) sign "+" and "— " correspond to the slow and fast MHD waves. It is worthwhile to note that Eqs. (10) and 
(11) are, of course, in agreement with previous results [1,4]. If we consider the isotropic pressure limit {p\\ = p± = p, 
C|| = c_L = c, and a = 5 = 0) and require = (? [4] then Eqs. (10) and (11) reduce to the form: 



1 %/5 %/5 

tan(/?(+0) = ±- and |m(+0)| = —^c = —^o, 



(12) 



as it was actually obtained in Ref.[4]. Whereas, applying the condition imposed in Ref.[l] a = c_l = \/2c|| we arrive 

at 



tan<^(+0) = ±^ and |u(+0)| = \/^c_l 

o V o 
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We must admit, here, that in the Eq.(9) for the modulus of the group velocity from the Ref.[l] factor of y | is missing. 

Apparently, substituting V{6 = 0) = ■ a = ■ ci_ and cos = \/l0/3 in the Eq.(9) (of the present paper) results 
in Eq.(13). This should not cause a confusion, as it does not affect the main result of the Ref. [1] — theoretical 
demonstration of possibility of the occurrence of the conical refraction in the anisotropic plasma plus calculation of 
the angle of the conical refraction, which is, actually, an experimentally observable physical quantity. We see that the 
results obtained in this paper are consistent with previous contributions [1,4]. However, it must emphasized that more 
general results given by Eqs.(lO) and (11) bring more diversion in the predicted phenomenon (see below). Besides, 
from technical point of view it should be much easier to produce the plasma satisfying condition (5) obtained in this 
paper, rather than condition a = c_l = \/2c|| requested in Ref.[l]. 



III. DISCUSSION 

The above consideration shows that when the wave normal N tends to reach a direction parallel to -Bo (i- c. when 
9 tends to zero), vector of group velocity u does not take the same direction, but its limiting position (when ^ — > 0) is 
the surface of an empty cone with opening angle 2(p{+Q) = 2 arctan[(/?(+0)] given by Eq.(lO). When 6 is zero exactly 
(which is quite unlikely from the experimental point of view as all waves produced in a laboratory may have only 
finite angle of aperture) then u is strictly parallel to Bq. The discontinuous behavior exhibited by the vector u can be 
explained by means of existence of so called [1,4] angular point which in the case of the anisotropic plasma revealed 
to exist when + — 2(q; + l)c| = 0. Thus, when + c^^ — 2(q; + l)c| = the angular point on the surfaces formed 
by the phase velocities (both fast and slow MHD waves) is located in the point of intersection of the ones with the 
direction of the magnetic field Bq. Which, in turn leads to the discontinuous behaviour of u and causes the effect. In 
other words, as it was mentioned in Ref.[l], the angular point exists when phase velocities of the slow and fast MHD 
waves do coincide. We observe form Eq. (6) [again, which was obtained from general dispersion relation (4) imposing 
condition (5)] that when 6 ^ Q slow ("— ") and fast ("+") branches of phase velocity equal to each other (viz., the 
square root vanishes) , which results in an uncertain behavior of its derivative by k i.e. group velocity. 

It is interesting to note that in differ to the result of Ref.[l], now imposing less restrictive condition (5) the angle of 
the conical refraction is not a fixed number but it depends [via Eq.(lO)] on the ratio ji?^. Which, in turn, depends 
on the plasma medium physical parameters such as density and the magnitude of the exerted external magnetic field 
Bo. 

As regards possible experimental corroboration of the effect, we may state the following: consider two collisionless 
plasma media contained in a miiform external magnetic field Bq . The media are separated by means of the insulating 
plane perpendicular to the Bq. Let us choose physical parameters in the second medium as to satisfy the condition 
(5). If we place a diaphragm between the first and second media then a narrow beam of MHD waves propagating 
from the first medium towards the second one, just after passing the diaphragm will deflect from the direction of Bq 
and form an empty cone with angle of opening 2(fi{+0). 

It is worthwhile to note that one of the possible applications of the effect of the conical refraction may be its usage 
in plasma diagnostics. The condition (5) explicitly depends on the magnitude of external magnetic field Bq. So, if in 
the above proposed experimental setup, in the second medium the condition (5) is not satisfied a detector of MHD 
waves placed just opposite to the diaphragm (see Fig.l) will fairly detect the waves. However, with variation of the 
magnitude of the exerted magnetic field we may achieve fulfillment of the condition (5) which will result in a steep 
drop of the detector's indication, as due to the effect of the conical refraction the MHD waves will deflect from the 
^ = axis and form an empty cone with opening angle 2ip{+0), thus avoiding arrival of the waves into the detector 
(see Fig. 2). Since the condition (5) explicitly depends on the external magnetic field, in the event of occurrence 
of the phenomenon it must be possible to measure its magnitude if the thermodynamical parameters (pressure and 
density) of the plasma are known or vice versa. Therefore, we conclude that presumably this effect can be used in 
plasma diagnostics. 

Besides, of a purely fundamental interest which experimental demonstration of the phenomenon might have, prob- 
ably, this effect will find its applications in practice, since through variation of external parameters (viz., varying the 
ratio a^/c^J it is possible to regulate opening angle of the empty conical beam of the MHD waves. Which means, 
that, for instance, with simple from technical point of view variation of the magnitude of the external magnetic field 
it is possible to regulate a direction of propagation of the energy carried by the MHD waves. Because, it is known 
that the energy propagates along the direction of the wave's group velocity vector. Therefore, one should not exclude 
possible applications of this novel effect (revealed by the general consideration employed in this paper) in the plasma 
devices like MHD generators are. 
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Finally, wc would like to emphasize on principal difference between the effect predicted in this paper and the effect 
of the conical refraction occuring in the optically and acoustically biaxial crystals. First, the most evident difference 
is that the phenomenon can take place in media with qualitatively different physical properties (magnetized plasma 
with isotropic [4] or anisotropic [1] pressure on the one hand and the biaxial crystals on the other). Second, in the case 
of biaxial crystals, let us take for concretness optical conical refraction, mathematical basis for theoretical formulation 
of the effect constitutes Maxwell's equations along with material equations and in addition electrical anisotropy of the 
medium is assumed, i.e. such substances (crystals) arc considered whose electrical excitations depend on the direction 
of the electric field [5] . Whereas in our case anisotropy in the medium is brought by the external magnetic field and 
Eulcr's equation for the magnetized, perfectly conducting plasma is used. The only similarity of these two phenomena 
is that the equations which relate the phase velocity of the waves with their direction of propagation in respect to the 
anisotropy axis (optical axis in one case and the external magenteic field in the other) — Presnel's equation of wave 
normals [5] and our dispcrtion relation (4) (including its special subcases [1] and [4]) have the same mathematical 
structure. Thus, geometrical shape of the surface formed by the group velocity vectors when the phenomenon takes 
place is qualitatively the same — surface of the empty cone. 
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MEDIUM I 




MEDIUM II 



FIG. 1. Sketch of experimental set up, which may 
be used for demonstration of the effect of the coni- 
cal refraction (see detailes in the text). The line with 
arrows represents the narrow, parallel beam of MHD 
waves. Note, that when in the medium II the condition 
a'^ +Cx — 2(a + l)c| =0 is not satisfied the beam goes into 
detector without experiencing any changes while propa- 
gating through the second medium. 



a^ + cl- 2{a + l)cl ^ 


+ c\- 2{a + l)c| = 


(2^(+0) DETECTOR 


MEDIUM I 


MEDIUM II 



FIG. 2. The same as in Fig.l. Note, that when in 
the medium II the condition + c\ — 2(a + l)cy = 
is fulfilled, the beam due to the effect of the 
conical refraction deflects from its initial direction 
(6 = 0) and forms an empty cone with opening angle 
2(p(+0) = 2 arctan[(^(+0)], thus avoiding arrival of the 
MHD waves into detector's receiver. 
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